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ABSTRACT 

(2) 

We define the notion of C^ +1 Ruijsenaars-Schneider models and construct their Lax 
formulation. They are obtained by a particular folding of the A 2 n+i systems. Their 
commuting Hamiltonians are linear combinations of Koornwinder-van Diejen "exter- 
nal fields" Ruijsenaars-Schneider models, for specific values of the exponential one- 
body couplings but with the most general 2 double-poles structure as opposed to the 
formerly studied BCjy case. Extensions to the elliptic potentials are briefly discussed. 
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1 Introduction 



We wish to present here an explicit Lax formulation for a subclass of one-body ex- 
tensions of the classical integrable Ruijsenaars-Schneider (RS) models quantized 
in 0. The quantum version of these extensions was formulated initially by van 
Diejen [§, || starting from the pioneering works of Koornwinder on algebraic for- 
mulations of extensions of MacDonald polynomials ||. This formulation, and the 
proof of quantum integrability relied upon analytical arguments using the newly con- 
structed Koornwinder polynomials as a basis for the wave function Hilbert space ||. 
A quantum formulation for the elliptic hamiltonians, conjectured in ||, was given in a 
series of works by Komori et al. H] using a corner-transfer matrix method combining 
functional representations of both bulk quantum i?-matrix and boundary reflection 
i^-matrices. General one-body extensions of the difference RS operators were thus 
built, limits of which could be identified with the one-body extensions of differential 
Calogero-Moser (CM) Hamiltonians constructed by Inozemtsev || and associated to 
BC N lattices [|. 

It must be noticed at this point that no clear classical limit of this construction 
exists, entailing as it does a coordinate-permutation operator in the functional rep- 
resentation of the quantum bulk i?-matrix and a coordinate-reflection operator in 
the quantum reflection K- matrix. In fact, no full classical Lax formulation exists for 
the most general Koornwinder-van Diejen (KvD) Hamiltonians; a first step in this 
direction was our identification of the BCn Ruijsenaars-Schneider Hamiltonians as 
linear combinations of particular KvD hamiltonians |I0| . In addition the classical 
r-matrix was obtained in this case, providing an interesting example of a dynamical 
dependance in both sets of dynamical variables, rapidities and positions. 

Our purpose is to describe a classical Lax formulation of a more general subclass 
of KvD potentials; more precisely the one-body part of the potential will exhibit 
the same 2 double-pole dependence as the full KvD; the coupling constants however, 
identified with the residues, still depend on one single coupling. This however marks 
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a qualitative advance compared to the previous construction of BCn RS model, since 
one there had only one double and one single pole. 

As in the previous case we will rely on a consistent Z2 folding procedure applied 
to a suitable A n RS model to obtain new classical integrable systems and their Lax 
formulation. The first suggestion of constructing (in that case) BCn and Cn RS 
models by Z2 folding of A2N(+i) RS model came from Ruijsenaars himself [jllj and 



was later explicited at the level of the Lax formulation |T!| . A related programme 
was applied to the simpler case of CM models 0, using folding procedures on the 
phase-space variables and twistings of the underlying Lie algebras to obtain general 
one-body extensions (see also 0). 

We here establish that this more general class of KvD potentials H (restricted 
to the hyperbolic case for the sake of simplicity) is obtained as linear combinations 
of pure RS Hamiltonians which may be described in terms of the root lattice of the 
twisted affine superalgebra OSp{2\2N) tw or C N+1 . The lattice characterizing the 
form of the potential indeed exhibits the shift of shortest roots by \ the derivation, 
characteristic of the Cffi+i root lattice |14| 0. At this time however we lack a deeper 



interpretation of the occurrence of this particular lattice. 

The problem of finding the Lax formulation for the most general four-coupling 
hyperbolic KvD potentials remains open at this point, but we conjecture that this 
subclass, exhibiting as it does the full pole structure, constitutes the best starting 
point for eventual achievement of this construction. 

The detailed analysis of elliptic generalizations will be left for later studies. How- 
ever it is established here that the most general 8 pole structure found in J7| (degen- 
erating to a 4 double-pole structure in the classical limit) may be similarly obtained 
by a generalized folding operation and the previous conjecture on the hyperbolic case 
is therefore extendable to the elliptic case. 

(2) 

A final comment regarding the denomination U C N+1 RS model" used here. The 



3 We are indebted to Luc Frappat for providing this reference and the identification of C^ +1 
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qualitative difference between the class of potentials obtained from the initial BCn 

(2) 

folding and this C^ +1 folding, i.e. the occurrence of a 2 double pole structure in 
the one-body term instead of a double-pole times single-pole structure, vanishes in 
the CM limit where poles add instead of multiply. It follows that the distinction 

(2) 

between BC N and C^+i CM models is essentially non-existent at least at the level 
of Hamiltonians, and the litterature (see for instance |§) rightly characterizes only 
"BC N type" models. 



2 The Cff +1 Ruij senaars-Schneider model 
2.1 Invariant subspaces of the A n RS dynamics 

(2) 

Presentation of the C N+l RS model is easily done in a more general framework where 
one constructs subsets of the phase space for an A n RS hierarchy so that they be 
invariant under a subclass of the flows generated by some specific combinations of the 
Hamiltonians. We now recall the definition of the A n RS hierarchy, for any n. 

The canonical variables are a set of rapidities {8i,i = 1 • • -n + 1} and conjugate 
positions qi such that {9i,qj} = dij. The Hamiltonians are initially defined as: 

h= J2 e ~ /361 ^ where fi = H <7<0 1/2 and 6 I = J2 6 *- 

IC{h-n+l},\I\=l i€l,k&I iel 

Function / may take different forms (rational, hyperbolic, trigonometric), resp: 

/(,) = !-?!; /(,) = l-ii^2; = 

q z smh -f sm 2 

The most general elliptic case where f(q) = (A + vV{q)) , V being Weierstrass func- 
tion will not be fully explicited here. Our construction may however be extended to 
it with due modifications to be discussed in the appropriate place (see Remark 2 in 
this section and Comment 2 in Section 5). 

The trigonometric and hyperbolic cases define the same model at least locally 
up to a redefinition of the parameters (the global structure of trigonometric vs. hy- 
perbolic models is however quite different, owing to qualitatively distinct topological 



properties, as can be seen for instance in [|15|). The rational case is obtained by an 
easy limit procedure from one of these models. We shall therefore consider in the 
following only the hyperbolic model. 

Let us note that f(q) = v(q)v(-q), with : v(q) = ^if^ = X' 1/2 jE^- These 
functions are thus rational fonctions of an exponential variable, z = e uq , defining 
A = e~ 27 . The square root is defined here with a cut on the real negative axis. 

Both functions / and v are periodic with period T = One introduces now a 
better adapted (although possibly degenerated) set of Hamiltonians as: 



K n = K+l + K+1 -iovl = l-- 



n + 1 



Ki = h h i^ = ^ = $>^// (1) 



K+i 



The negative- index Hamiltonians are in fact obtained from traces of negative 
powers of the RS Lax matrix, using its remarkable Cauchy structure [|TJ to con- 
nect L~ l with I}. One then considers any idempotent bijection o over the index 
set {l---n + 1}, separating them into invariant singlets o~(i) = i and doublets 
a{i) = j , a(j) = i , j ^ i. It is easily shown that: 

Proposition 1 For any idempotent bijection a, the manifolds defined by: 

Vi, qi + q^t) = 0(T) ; 9i + 9 a{i) = 0(— ) (2) 
are kept invariant by the evolution generated by the Hamiltonians Kl+ 2 ~ l = ■ 

e Wi _ e -l39i 

Proof: The evolution equations are: {qj,Tti} = (3 / , ^ // an d {^j^l} = 

E 2 fl E^/ 1/2 )'(^-^)-E 2 h E^/ 1/2 )'(^-^)- o) 

Invariance of the manifolds (0) under (|3|) is straightforwardly obtained from the 
parity and T-periodicity of the two-body potential function /; the explicit 9 parity and 
^-periodicity of Hi; and the bijectivity of a allowing adequate index redefinitions. 

The number of a-stable indices is a priori arbitrary. However if there are more 
than two such indices, one will unavoidably have exact equality mod. T of at least two 



position variables q and the invariant submanifold will actually lie in the singularity 
hyperplanes qi = qj(T). This case must therefore be eliminated from our discussion, 
which leaves us with only three possibilities. 

Case 1: For n even (odd number of sites), one may only consider the case of one 
stable index. This construction gives us the BCtl RS model |]TT , 12] , [10]. 

Case 2: For n odd (even number of sites), one may first consider the case of zero 
stable index. This leads to the Cn+i case [[II], [L3|. 

Case 3: For n odd (even number of sites), one may then consider the case of two 
stable indices hereafter denoted and 0. The only non-singular choice for values of the 
two fixed position coordinates is then q = and gg = T/2 up to trivial permutation. 

(2) 

Connection to the superalgebra C N+1 with n = 2N + 1 will be established presently. 

Remark 1: Regarding the rapidities, we will here restrict ourselves for the sake 
of simplicity to the choice 9 = 6*o = 0(^p). This implies immediately that the 
restriction of 7i to @ is identically 1. We shall comment at the end (Comment 3, 
Section 5) on the implications of other possible choices. 

We now establish the important properties of this consistent restriction. We 
first introduce a better adapted set of conjugate variables (for the full phase space) 
defined to be {^(g* - q^)),^(0i ~ {75 (ft + 9<r(<)) = 9o 73^ + ^«) = 

, V« , a(i) ^ i; finally q = qfr, 9 = 6fr, g 5 - f = q$, 0q = 0jj"- 
Corollary la: These submanifolds are obviously endowed with a symplectic form 
{ }rest where conjugate variables are respectively the restrictions of 75 (ft — Qa{i)) an d 
4^(^i — for i such that a(i) 7^ i. 

Corollary lb: The Hamiltonians Tii, restricted to the submanifolds (||) endowed 
with the symplectic form { } rost build an integrable hierarchy. 

Proof: An obvious rewriting of Proposition 1 states that the Hamiltonians Tii ex- 
panded around the values Of = 0, qf = 0, have no linear term in Of, qf: indeed such 
terms would trigger a non-trivial dynamics of 0f,qf on (U). This leads us to define 
U\ est = " ftx(i)), ^=(Pi - 6*®W = OX = 0) and now H t = 71^+ 

quadratic terms in qf-, Of 
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In addition the Tii Poisson-commute by construction. This is now expressed as : 
{Ha, ?4}fuii = = {H. r a est , Hl est } lu i\ + 6^ x h qf x • • • 

which implies in particular: {H r a est , H[ est } {u u = 0. 

Finally one has {H r a es \ H£ est } m = {H r a est , Hl est } vcst = by definition of the re- 
stricted symplectic structure on the submanifolds, thereby ending the proof. 

For the sake of simplicity we shall from now on drop the "rest" index in HCi. 

Remark 2: In the elliptic case the potential function is biperiodic, hence one may 
allow for 4 fixed coordinates since two independent periods are available to define the 
invariant subspaces (replacing the single period congruence parameter T in @). 

We now formulate: 

Proposition 2 The Hamiltonians K\ and K_\ are identical, and therefore equal to 
Hi, on the submanifolds ([Jj. 

The proof is a trivial consequence of the bijectivity of o~\ the fact that on the 
reduced submanifolds (H) one has up to a full period Q a ^\ = —6i and q a ^ = —qi for 
all indices; and the parity and periodicity properties of the potential /. 

The immediate crucial consequence is that the restriction of the Hamiltonians Tii 
can be rewritten as: Hi = Pi(TrL m ) where P; is the l-th Newton polynomial (sum 
of rank I minors) and L is the restriction of the Lax matrix of A2N+1 RS model to 
the submanifolds (0). We are therefore provided with a Lax representation for the 
integrable hierarchy "Hi. 

We now discuss the particular dependance of the potential terms in the Hamil- 
tonians, in order to justify the claimed connection to the twisted affine superalgebra 

(2) 

C]v+i • Specific connection is the following: the fixing of gg to the value T/2 introduces 

an imaginary shift in the lattice describing the position dependance of the potential 

(2) 

function, turning it into the root lattice of C N+1 . 

More precisely the original dependance of the potential function was in the vari- 
ables qi — qj, i,j G {1 ■ • ■ 2N + 2}, associating it with the root lattice of A 2 n+i- After 
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this particular folding we now get a dependence in qi — qj)qi + qj;qi,2qi;qi + T/2;i,j G 
{1 • • • iV}. This immediately leads to seek for a lattice with three root lengths together 
with the simultaneous existence of shifted and unshifted shortest roots. By explo- 
ration |14[ the only possibility is C^' +1 . 

These variables may now be rexpressed as scalar products < a.Q > where Q is 
the (N + l)-component vector in the dual of the Cartan algebra with coordinates 
Qi ' ' ' Qn on the dual of the Cartan algebra of the simple Lie algebra and T along the 
direction of the derivation generator d of the full Cartan algebra; a is any root of 
Cjp +1 [14 1, exhibiting in particular the shift of the shortest roots by ~ the derivation 
generator characteristic of C N+1 . It remains at this time still a formal connection, 
and we have no interpretation of the occurence of a twisted affine superalgebra in this 
context. 

The pure BCn folding (case 1) differs qualitatively in that the one-body part of the 
potential for the first Hamiltonian (i.e. linearly dependent upon single exponentials 
of rapidities) does not contain a double pole at half period due to the absence of the 
extra dependance in qi + T/2. In the non-relativistic CM limit however, the one- 
body potential in the first (quadratic) Hamiltonian obtained by both foldings exhibits 
double poles at half-period (occurring from the folding of terms sinh(gj — qj)^ 2 at 
qj = —qi giving in particular a term cosh(gj)~ 2 by doubling of the sinh; or at qj = ^); 
and double poles at integer period (occurring from the folding of terms sinh(gj — qj)~ 2 
at qj = 0). The technical reason is the multiplicative nature of the potential terms 
in a RS model as opposed to the additive nature thereof in a CM model, and is thus 
related to the fundamental difference between RS models realized [16| on a Heisenberg 



double |17|, [19[ and CM models realized on a cotangent bundle ||20|| . Actually, the 

(2) 

only difference between BCn and C N+1 foldings of CM models lies in the value of the 
residues (coupling constants) which in any case may be extended to take any value, 
hence it turns out to be irrelevant. 
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(2) 

3 The C N+1 Lax operator and r-matrix 

We shall from now on use the notation — % instead of cr(i), defining that —0 = 0. The 
Lax formulation of C^ +l RS system is therefore obtained as a folding of the A 2 tv + i 
TIJI explicitely built as follow: 



case 



We first label the 2N + 2 rapidities i = —N • • • 0, 0, • • • N} and conjugate posi- 
tions {qi, i = —N • • • 0, 0, • • • N}. Independent phase space variables on the restricted 
manifolds are here chosen to be qi,6i,i G {1---N}. Note that on the restricted 
manifold the Poisson structure {} res t in terms of these variables reads qj} = ^Sij 
which will be responsible for an overall | factor in the Poisson structure. This fac- 
tor will be omitted from now on, corresponding to a normalization of the Poisson 
bracket as {Oi,qj} = <^j. One thus identifies 6{ = £iO\i\ and qi = £%q\i\ + #i,of with 
Vi G {1 • • -N}, Si = 1 = — E-i and Eq = Eq = 0. 

N 

The Lax matrix for the A 2 n+i cases reads: £ = £jj where Cij(q-N, • <?v? 0j) 

i,j=-N 

c(qi — qj) e~P dj fyy, {e^} is the standard basis for (2N + 2) x (2iV + 2) matrices and 
^) = id^ = (l-A)S- 

The Lax matrix for the C$ +1 Ruijsenaars-Schneider systems then reads: 

N rp 

L= ^ L ij e ij With L ij = £ij(-QN, ■■; "91,0, — , ?1, q N , Ej 0\j\) (4) 
i,j=-N 

It is now well-known that the Lax operator £ satisfies the quadratic fundamental 



Poisson bracket f2~l|| : {C f £} = C (g> C a x — a 2 C <E> C + C 2 Si L\ — L\ s 2 C 2 , where 



L\ = C®1, C 2 = 1®C and the quadratic structure coefficients read: a\ = a + w, Si = 

N 

s — w, a 2 = a + s — s n — w, s 2 = s^+w. As usual, for any matrix M = M^ki &ij®e k i 

ijkl=-N 

the operation n is defined by: M[j kl = M k uj. Matrices a, s, w take the form : 

N N N 

a = -«2J a jk ejk ®e k j, s = s jk e jk <S> e k k, w = a jk ejj <g> e fcA: with at = /3 - 

j,k=—N j,k=-N j 7 k=~N 

and Oj fe = (1 - dj, fe ) coth~(g 3 - - q k ) = (1 - d ijfe ) , s jfe - 



2 "" J ' Zj-Zk wah%(qj — q k ) 



Remember that the most general structure of Poisson bracket for a Lax operator 
of a Liouville-integrable system is a linear one [ffiql : {L f L} = [r, Li] — [r^, L2]. The 
quadratic form corresponds to the general case |23|] where the r-matrix itself assumes 
a linear dependence in L of form: r = b L 2 + L 2 c with 6 and c arbitrary matrices 
yielding : a\ = c n — c, a 2 = b n — b, Si = c + 6 71 " and s 2 = . 

We now straightforwardly extend the previous computation of the BCn 



Ruijsenaars- Schneider r-matrix. The Lax operator also satisfies a quadratic fun- 
damental Poisson bracket, again exhibiting the dependence of the structure matrices 
a and s on both sets of dynamical variables. The occurrence of two zero-type in- 
dices and requires extra caution at one particular place, eventually leading to 
supplementary signs in the formula for factors of the r-matrix. Note finally that this 
derivation is equivalently applicable to the Cjv-type RS Lax matrix, this time by 
altogether removing the zero-indices. 

The r-matrix structure is again completely defined by a quadratic Poisson bracket 
with a,i, a 2 , Si and s 2 changed into: a x — > 5i = ax, a 2 — > a 2 = 5i + s\ — s 2 = 
a 2 + T 7r - r, si — > Si = Si + a + , s 2 — >■ s 2 = s\ = s 2 + a w + r = s 2 + a + r, with: 

^ it ^ .4 

1" = «^] £I« (-^~ £ fc-i) e«®e w , a = a ^ (5^ Sj - (1 - ^j) a i3 -) --p e^e^-i, 

i,k,l=-N ' i,j=-N % 

1J2 . , ^ t v 1 \,\ ^i + A^j 1 + A t m j + tj m 
defining = coth(-(% - g,) + 7) = _ r , s< = —7 + }^ 9 ' Ui = 



m=-N 



f> ifc + t w - and A=f\ § XZk )^ < " '"■< ' 

k=-N k = -N, 

k i 

Occurrence of the final sign factor e i7r<5i >° is the only qualitative modification in- 
duced by the supplementary fixed index 0. In addition the full r-matrix gets an 
overall | factor due to the normalization of the restricted Poisson bracket. 

As in the BCn case this quadratic r- matrix structure is fully dynamical, depend- 
ing both on the positions q^s and rapidities Ots. 
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4 Koornwinder-van Diejen versus Cjp +1 Ruijsenaars- 
Schneider Hamiltonians 

As in the BCn case the Hamiltonians Hi generated by traces of powers of the 
Cjy +1 Lax matrix (|j) can be rexpressed in a more interesting form as: 



Hi = u J c ,i-\J\ e ~ P6eJ Uj, with eJ = {sj\j\,j G J} and 

j c {1..JV}, |j| < ; 

ej = ±1, j e J 

w ** e n / i/2 ^-^)= e n ( 5 ) 

s c Ar- = K U --8T U{o>5} s es sc.Ak sss 

S = -S, |S| =p fee Ak\S S = -S, \S\=p k e ^4k\S 

We now recall the form of classical Koornwinder-van Diejen Hamiltonians ||: 



E ^-ui V *J*J° V %°> ( 6 ) 

J C {l..iV}, J| < i 
Ej = ±1, JEJ 



where, after some rearrangements required to reintroduce the indices 0, 0: 



VeJ > K = Y[ Wr ( q i"> EI v(qj -q k ), with. w r (qj) 



jeeJ , v{2q i )v{q j )v{q j - f) 



fe S -4 A' U - £ J 



and u K , P = (-iy Yi n^(») n ! Vr-^r n w 

eJ C Ak l ^ £l i,i r £ el i € el 



The potentials w are particular functions explicited in and may be interpreted 
physically as an interaction with some external field. 

Direct computation yields: V £ j-jc V- £ j-jc = J^J w r (qj) w r (—qj) f^j. 

jeeJ 

Setting w r (qj) = 1, that is w(qj) = v (2 qj) v(qj) v(qj — which is an admissible 
choice according to |J, Hi (|5]) takes actually the same form as Hi (^), up to the 
crucial change of Uk,p into Uk, p - They are generally not equal, except for p = 0, 
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where trivially: Uk,o = 1 = Uk,o- When p = 1, one gets: 

-i— r -i— r J 1 

i6^ K \{0,0} ke 

fc^i fc^O fc # 



Evaluation of a suitable contour integral as in gives the Liouville-type functional 

■7(2|g| 
sinh7 



identity: £ J] - q k ) = rinh ^ +2 > , that is: ^ = - 



We now recall the general theorem established in [|l(| : 



Theorem 1 Let and 9{, i 6 N, 6e a sei o/ conjugated variables such that {6t, qj} = 
5ij. Let I and K be arbitrary finite sets of indices included in N. Assume the existence 
of a set of complex functions Uk, p depending upon the set of indices K and a natural 
integer p, and of another set of complex functions v £ jj depending upon the sets of 
indices J and I (J C I) and a \ J\-uple of signs e = (ej, j G J), such that: 

• uk, p and v £ jj be independent of the rapidities 9iS. 

• u K ,o = 1, ^0,7 = 1, and v^yj ^ 0. 

• S 1 = {h\ = ujc^_\j\ e~ /36eJ v £ jj, I G {l..|/|} } be a family of Poisson- 

J C I, | J\ < l 

s 3 = ±1, 3 e J 

commuting functions (9 e j = 

If there exists a second set of complex functions uk, p obeying the first two condi- 
tions; such that S 1 = { h\ = mj c ,z-|j| e~ f39eJ v e jj, I G {l..|/|} } be a new 

J C I, \J\ < l 

Ej = ±1, j S J 

family of Pois son- commuting functions; and uk,i = uk,i + c i(|-^|); then there exist 
coefficients c r (m), (r, m) G N 2 , independent of all dynamical variables, connecting the 
two families of Hamiltonians as: hj = X^=o c i-s(\I\ — s) h r s , with Vm G N, co{m) = 1. 

Hence the two relations for p = and p = 1 between the Uk, p 's and Uk^s are 
sufficient to establish that the two sets of Hamiltonians define the same family of 
commuting dynamical flows, namely one set of Hamiltonians is a triangular linear 
combination of the other set. 



11 



5 Comments 

(2) 

1. It follows from the uniqueness theorem that the C N+l RS hierarchy defined here is 
equivalent to the KvD hierarchy for a particular set of couplings /i, /i, ^, ^. However 
as already emphasized the full KvD pole structure is now obtained, contrary to the 
pure BCn case where one coupling constant is actually set to zero. To obtain the 
complete KvD set of hyperbolic potentials one clearly needs to define an "extension" 
(with the same meaning as the "extension" of BCn CM models in |9], |J leading to 
the full Inozemtsev potentials) of the C N+1 Lax formulation. We hope to report on 
this problem soon. 

2. In the elliptic case where 4 coordinates q are fixed to various combinations of the 
two half-periods, the folding leads to a one-body potential with 4 double poles located 
at qi— (half-integer linear combinations of the two periods): indeed the term depend- 
ing on 2qi leads to a duplication of the 4 poles induced by the 4 fixed coordinates. 

This is the correct denumbering and location of poles in the classical limit of 
the quantum general elliptic Hamiltonian constructed in 0, where 4 pairs of poles 
separated by an order h degenerate to these 4 double poles. Of course, here again 
one does not obtain the full set of coupling constants (or equivalently the residues 
at the poles). The conjecture that folded RS models are the correct starting points 
for construction of a full KvD-Hikami-Komori classical Lax-type representation may 
therefore be extended to the elliptic case, with a suitably generalized underlying root 
lattice. More precisely, we conjecture that this new root lattice is associated with 
twisted toroidal superalgebras. Indeed it exhibits three root lengths (pointing to 
superalgebras) with half-period shifts (pointing to twisted algebras), precisely two 
independent shifts of the shortest roots, interpreted as corresponding to the two 
independent derivations characteristic of toroidal algebras. 

3. We have restricted ourselves here to the choice 6q = 8q = 0. In fact one can also 
choose either or both to be equal to a half-period jj-. The corresponding phase space 
manifold is also invariant under the Hamitonians Tii. These Hamiltonians therefore 
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define an integrable hierarchy. At first sight it is not identical to the one defined here: 

the Hamiltonians indeed exhibit modifications of several relative signs in the potential 

terms; for instance the pure potential in the first Hamiltonian will become Uk,i = 

T 

J { v (Qk) ~ | J v (-^ ~Qk) if #o — lJ k- Whether these hierarchies are genuinely new 

fe 6 Ak k e A K 

k k ^ 

(2) 

integrable systems or may be obtained from our original C N+1 Hamiltonians by some 
redefinition of parameters is an open question which we will not address here. 
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